Scattering of a Dirac electron on a mass barrier 
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The interaction of a wave packet (and in particular the wave front) with a mass barrier is in- 
vestigated in one dimension. We discuss the main features of the wave packet that are inherent to 
two-dimensional wave packets, such as compression during reflection, penetration in the case when 
the energy is lower than the height of the barrier, waving tails, precursors and the retardation of the 
reflected and penetrated wave packets. These features depend on the wave packet envelope function 
which we demonstrate by considering the case of a rectangular wave packet with sharp front and 
trailing edges and a smooth Gaussian wave packet. The method of Fourier integral for obtaining 
the non-stationary solutions is used. 
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I. INTRODUCTION 

Progress in nanotechnology has triggered a broad in- 
terest in low dimensional physics. The separation of 
graphene by mechanical exfoliation in 2004i^ opened a 
broad range of activity for theoretical and experimental 
researchers. The interest rose due to the peculiar prop- 
erties of graphene such as ultra-relativistic behavior of 
charge carriers with a Fermi velocity 300 times smaller 
than the velocity of light, the linear spectrum close to 
K and K' points in the momentum plane that can be 
described by the massless Dirac- Weyl equation^, an un- 
conventional quantum Hall effect^, and the perfect trans- 
mission through arbitrarily high and wide barriers, the 
so called Klein tunneling^. 

The energy levels or the spectrum of the electron sys- 
tem typically investigated as being the main character- 
istics of quantum nano-structures (say, such as quantum 
dots) are found by considering the stationary Schrodinger 
equation^. Recently, a scientific interest has shifted to 
the investigation of quantum dynamics including quasi- 
bound states^, electron beams, wave packets, and their 
control by means of barriers and other nonhomogeneous 
structures^—. Unfortunately due to above mentioned 
Klein effect and the gapless spectrum the control of elec- 
trons in graphene by means of electric fields is inefficient. 
Therefore, the search for other possibilities to control 
electrons in graphene become of interest. Recently, the 
creation of a gap in the electron spectrum^ as well as 
the control of the valley isospi n 12 i 13 was predicted by in- 
troducing a mass term into the Dirac- Weyl Hamiltonian. 
This possibility was experimentally demonstrated by the 
proper arrangement of dopants in the graphene sheets 
or by inducing electron-electron interactions^. There- 
fore, the investigation of the non-stationary solutions of 
the Dirac equation describing the interaction of electrons 
with mass barriers becomes timely. Such sophisticated 
systems and problems are a challenge for numerical sim- 
ulation. But in order to obtain a physical understanding 
of the behavior of such systems it is very helpful to invent 



and analyze simple models that are able to demonstrate 
the main physical features of interest. 



The purpose of the present work is to demonstrate the 
mean properties of wave packets, such as their reflection 
and penetration into mass steps and barriers, the appear- 
ance of precursors, the formation of evanescent waves 
and the zitterbewegung, using very simple analytically 
solvable one-dimensional (ID) models. Comparing the 
results obtained for a rectangular wave packet with the 
ones having a rather soft Gaussian shaped wave packet 
we demonstrate that their propagation depends on the 
shape of the wave packet. In our analysis we pay special 
attention to the appearance of so called precursors that 
were predicted by Sommerfeld nearly a hundred years 
ago for electromagnetic waves travelling in a dispersive 
medial. We also focus on the waving after-effects that 
can be related to the zitterbewegung that was predicted 
by Schrodinger for relativistic electrons^, and which we 
predict should be observable in graphene^. The results 
were obtained applying the Fourier integral technique for 
solving the non-stationary Dirac equation, and analyzing 
the solution by means of integration in the complex wave 
vector plane. 



The paper is organized as follows. In Sec. II we in- 
troduce the problem and the method of its solution con- 
sidering the reflection of the Dirac wave packet from a 
hard wall. In Sec. Ill results for scattering of a plane 
wave type wave function by the mass barrier is discussed 
which enables us to present the non-stationary solutions 
for the reflected front in the form of a complex integral in 
Sec. IV. The penetration of the front of the wave packet 
into the barrier is considered in Sec. V, and in Sec. VI 
the scattering of the rectangular wave packet is studied. 
Sec. VII is devoted to the description of the scattering of 
a Gaussian wave packet. Our conclusions are presented 
in Sec. VIII. 
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II. REFLECTION OF A WAVE PACKET FROM 
A HARD WALL 

We consider the motion of a ID wave packet that is 
described by the following Dirac-Weyl equation: 



d 



d 



-ia x —^{x,t) + Q{x)a z ^{x,t), (1) 

where a x and <j z are the Pauli matrices, and the symbol 
0(x) stands for the Heaviside step function whose value 
is zero for negative argument and unity for positive one. 
The last term will be referred to as the mass term, that 
leads to a gap in the electron spectrum close to the K 
and K' points. 

In order to simplify the notations the above equation, 
the results are presented in dimensionless notations based 
on the height of the mass barrier V. So, the energy is 
measured in V units, the time — in h/V units, and the 
unit of length is Hvf/V, where the symbol vp stands for 
the Fermi velocity. For the sake of illustration we took 
V = 53 meV, what is achieved by depositin g th e single 
graphene layer on a BN substrate (see Ref. Il9h . Then 
the time unit is a hundredth of ps, and the unit of length 
is about 10 nm. 

Presenting the wave function as a two component 
spinor 



V(x,t) 



u(x, t) 
v(x, t) 



(2) 



we have to solve the following set of two differential equa- 
tions for the wave function components: 



u t = -v x - i@(x)u, 
v t = —u x + iO(x)v. 



(3a) 
(3b) 



Apart of the wave function components themselves the 
wave packet can be characterized by the local density 



(4) 



p{x,t) = \u{x,t)\ 2 + \v(x,t)\ 2 , 

the local current 

j(x,t) =2Re[u*(x,t)v(x,t)], (5) 

and some averaged values: the number of particles in the 
wave packet (the norm of the wave function) 

/oo 
dxp(x,t), (6) 
-OO 

the mean position 

/oo 
dxxp(x,t), (7) 
-OO 

the mean velocity 

V(t)=N- 1 (t) dxj(x,t), (8) 



and the width of the wave packet 



(t) = N-\t) dx[x- X(t)} 2 p(x,t). (9) 



In order to introduce the necessary definitions and il- 
lustrate the method of solving the time dependent Dirac 
equation we start with the most simple problem: the re- 
flection of the wave packet from a hard wall. This means 
that Eqs. without the last term (i. e. the massless 
Dirac equation) will be solved in the negative part of the 
x-axis (—00 < x ^ 0) and the hard wall will be taken 
into account by the boundary condition 



u(0,t) + iv(0,t) = 



(10) 



which was derived and discussed in Ref. Il3l 

This solution can be easily composed of two (incident 
and reflected) freely propagating wave packets. The in- 
cident wave packet moving to the right side can be pre- 
sented as 

(x,t) = ^ J c lq{x ^^>(t-x), (11) 

that satisfies the above massless equation with any en- 
velope function $(x). In order to satisfy the boundary 
condition (fTU|) at any time t we have to add the reflected 
wave packet 

¥ lfi \x,t) = -J ^ 1 i ^ e- iq{x+t) t>(t + x), (12) 

that propagates to the left. 

The most interesting feature of the wave function con- 
structed in such a way 

<&(x,t) = ¥ in \x,t) + ¥ lfi \x,t) 1 (13) 

and describing the Dirac electron reflection from the hard 
wall, is the absence of interference of the incident and 
reflected wave packets. Indeed, denoting 

yr ± = e M±x-t)$( tTx ^ ( 14 ) 

we present the density as 

p( 2 ;,t) = |/ + -z/_| 2 + |/ + + l /_| 2 

= l/ + | 2 + l/-| 2 = l<!Ki-z)| 2 + l^ + z)| 2 , 
and the current as 

j(x,t) = He + )(/+ + «/-)} 

= Re {|/+| 2 -|r| 2 +2fiRc (/;/_)} (16) 

= \f + \ 2 -\.r\ 2 = Mt-x)\ 2 -Mt + x)\ 2 . 

Thus, the density and current are expressed just through 
the corresponding individual values of the incident and 
reflected wave packets without any interference terms. 
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In order to trace how the reflection depends on the 
form of the envelope function we consider in detail two 
extreme shapes of wave packets: a wave packet with a 
rectangular envelope as an example of a wave packet with 
abrupt edges, and a rather soft Gaussian one. 

We choose the following envelope for the rectangular 
wave packet: 

$(x) = -0(a 2 /4-x 2 ). (17) 
a 

Inserting this into (JTHH]) we obtain the following mean 
values 
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V(t) 
w(t) 
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(18a) 
(18b) 
(18c) 



characterizing the reflection of the rectangular wave 
packet from the hard wall. They are shown in Fig. [lja). 
We would like to draw your attention to the following fea- 
tures of this simple example of reflection: (1) the wave 
packet is reflected before it reaches the wall, (2) it is com- 
pressed during the reflection, and (3) the wave packet di- 
mensions are restored after the reflection. This behavior 
can be compared with the one of a rubber ball hitting a 
hard wall: due to its finite extension the ball changes the 
direction of motion before its center reaches the wall and 
it is squeezed during the collision. 

In the case of a Gaussian wave packet the following 
envelope function was chosen: 



$(x) 



(19) 



where the parameter a characterizes the width of the 
wave packet and plays the same role as the parameter a 
in the previous case of a rectangular wave packet. Now 
inserting this expression into Eqs. pTHT3|) and then into 
Eqs. (JTHH]) we obtain the following averaged values: 



(20a) 
(20b) 

(20c) 



X(t) 
V(t) 

W(t) = CT + t 2 - 



c -t /a* _ t . er f(i/V2?), 
-eri(t/V2o), 

l^ e -t 2 /2, _ t . erf (t/ ^ 



n 2 



where the symbol "erf" stands for the error function. 
These dependencies are shown in figure Fig. D^b). Com- 
paring both Figs, (a) and (b) we see that qualitatively 
they are the same, and consequently, the reflection of the 
Dirac wave packet from the infinite wall is not sensitive 
to the form of the envelope function. This non-sensitivity 
is quite expected because the dispersionless propagation 
of the wave packet described by Eq. ([lj without the mass 
term is not spoiled by the dispersionless boundary con- 
dition ([ID]). 
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Figure 1: (Color online) Evolution of mean values for the 
Dirac wave packet in time: the mean coordinate of the wave 
packet - green dashed curve, the mean velocity - blue dotted 
curve and the width - red solid curve, (a) - the rectangular 
wave packet with envelope (|17[) and (b) - the Gaussian one 
with envelope (I19|) . 



III. MASS BARRIER OF A FINITE HEIGHT 

Now we consider our main problem: the interaction of 
a wave packet with a mass barrier of finite height. This 
problem is much more complicated than the previous one 
because now the wave packet spends some time inside 
the barrier region which is a dispersive medium, and as a 
consequence the wave packet will not conserve its shape. 

As the Dirac equation is linear with coordinate inde- 
pendent coefficients a natural way to solve the problem 
is by using Fourier transformation. Thus, we choose the 
initial condition (say the position of the wave packet at 
time t = 0), expand it into Fourier series (here integral 
over momentum fc), and then change the exponent in 
the integrand by the function that obeys the time de- 
pendent Dirac equation in both regions (the barrier and 
free motion region) and is consistent with the boundary 
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conditions at the point x — 0. For this purpose we re- 
cap briefly the main results for scattering of a plane- wave 
(exponent) type wave function. 

We assume that in the region of free motion (—00 < 
x < 0) there are incident and reflected waves with energy 
e = k, and here the total wave function reads 

* (froc) (;M) = 1- (in) 0M) + # (lfl) (z,t) 




In the barrier region (0 < x < 00) there is only an outgo- 
ing wave moving to the right with the same energy k and 
momentum k = \Jk 2 — 1. We present its wave function 
as 



V {tI \x,t) = 



T 



VkTT 
2k \ y/k-1 



-ikt 



(22) 



Now equating both wave function components at the 
point x = we obtain the following wave reflection and 
penetration amplitudes: 

R{k) = k-K, (23a) 
T(k) = Vk (Vk + T - Vh~tj . (23b) 

Eqs. ((2TJ and j22]) together with definitions ([23} en- 
able us to present the time dependent wave function of 
any wave packet in the form of the integral along some 
contour in the complex fc-plane. 



IV. REFLECTION OF A STEEP FRONT 

First we consider the motion of the wave packet with 
rectangular envelope function. Fortunately, due to the 
linearity of the Dirac equation this problem can be de- 
composed into the motion of two fronts corresponding to 
the leading and trailing edges of this wave packet. That is 
why we start with the reflection of a steep front from the 
mass barrier, and choose the initial incident wave packet 
as 

f (inc) (a;,0) = -j= ( j e lqx O(-x). (24) 

It coincides with the freely propagating front where the 
leading edge has reached the barrier. This function can 
be presented by the following Fourier integral: 



where 



dxe [t{q ~ k)+a]x 



1 



i(q — k) + a 



(26) 



Here the symbol a stands for the regularization param- 
eter — a small positive value that will be set to zero at 
the end of the calculation. 

Now according to our strategy we replace the exponent 
exp(ifcx) in the integrand of the Fourier integral (|25|) by 
the solution of the time dependent Dirac equation (|21I22[) 
describing the reflection and penetration into the barrier 
of this exponent type wave function. The wave function 
of the wave packet obtained in this way consists of three 
parts. Two of them are defined in the free motion region 
(—00 < x < 0) and describe the incident front 



¥ in \x,t) = - 



1 

71 



1 

2ni 



dke ik(x-t) 
k — q + ia 



(27) 



and the reflected one 
^ rfl )(z,£) 



_L I 1 ^ _L f e~ lk{x+t) (k- K )dk (28) 



V2 \ -1 



2ni 



The third part of the wave function is defined in the 
region < x < 00 and describes the front that penetrates 
the barrier: 



^^\x,t) 
1 



2 3/: 



7" J C 



e t(KX-kt) dk 




(29) 



Integral (|2"T)) is trivial because its integrand has just 
a single singular point in the lower part of the complex 
fc-plane, namely, the pole at the point k = q — ia. Thus, 
choosing the integration contour C by passing that point 
from above and enclosing it by the upper or lower semi- 
circle (depending on the sign of the parameter x — t in 
the exponent) we calculate the residue at this pole and 
obtain the following wave function of the incident front: 



* (in )(x,i) 




0(t - x)c lq ( x -^ (30) 



that describes the motion of this front with a constant 
velocity equal to unity due to the absence of dispersion 
in the free motion region. 

The two other integrals (|2"5)) and (f2"9")l are more com- 
plicated because of the radicals k = y/k 2 — 1 in their 
integrands leading to branching points at k± = ±1 in 
the complex fc-plane. Consequently, we have to make a 
cut in this plane and choose the position of the integra- 
tion contour accordingly. How this is done by taking the 
causality principle into account is shown in Fig.[2ja). The 
choice of the Riemann sheet is defined by the requirement 
Im k > in the upper half-plane, and this is determined 
by the fact that at large momentum k (positive or neg- 
ative) the value of k should approximately coincide with 
the value of fc. The choice of the integration contour C 
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Figure 2: (Color online) The complex fc-plane, pole, cut and 
contour for the integration in equations (|28[) and (|29[) : (a) - 
initial contour, (b) - transformation of the contour in the case 
q > 1 and (c) - in the case q < 1. 



laying above both singularities (pole and cut) is in agree- 
ment with the requirement that the integral (|28|) should 
be zero in the case of x < — t as the reflected front can 
not move faster than with unit velocity. 

In the case of x > — t the argument of the exponent 
has opposite sign and that is why the contour can be 
enclosed by the lower semi-circle and transformed into 
two contours encircling the singularities as it is shown 
in Fig. [2lb,c), corresponding to the contributions of pole 
and cut. 

The pole contribution is calculated by means of the 
residue technique and it gives 



e(x + t)R(q)e- iq ( x+t) 



1 

71 



i 



-l 



(31) 



It coincides with the result presented in section Mil for 
scattering of a plane wave. The single difference is that 
now the exponent has a steep leading edge. 

The calculation of the cut contribution is more compli- 
cated. It can not be calculated analytically, and a numer- 
ical evaluation of the integral is necessary. This integral 
depends essentially on the energy e = q of the incident 
front, namely, whether it is larger or smaller than the 
height of the barrier. 



A. Above the barrier reflection (g > 1) 

For q > 1 the contours Co and Cb are separated hori- 
zontally (see in Fig. [DJb)) and the pole of the integrand 
doesn't complicate the calculation of the cut contribution 
along the contour Cb- That is why we transformed this 
contribution into two real integrals that depend on the 
single argument £ = x + t: 



EE 

— ^cut 



II 



k 2 



{— ycos(fc£) +isin(fc£;)} 



(32) 



and integrated this numerically. 

A typical result is shown in Fig. |3(a). There are two 
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Figure 3: (Color online) The wave function of the reflected 
front as a function of the variable £ = x + t: blue dashed curve 

- the real part of the component it = — v, green dotted curve 

- the imaginary part of it, red thick solid curve - the local 
density, and black thin solid curve - the same local density 
with only the pole contribution included; (a) - above the bar- 
rier reflection with energy q — 2 and (b) - below the barrier 
reflection with q = 0.6. 



points worth to be mentioned. First, at the beginning 



(i 



of the front there is some overshoot of the wave function 
modulus squared (the local density) as compared with 
the one when only the pole contribution is taken into ac- 
count (as was already mentioned the pole contribution 
coincides with the result of plane wave scattering). Sec- 
ond, the above mentioned local density tends to the pole 
contribution with increasing £ = x + t (when we are going 
away from the leading edge of the front) but notice that 
some oscillations are still present. Following^ the over- 
shoot can be explained as follows. According to Fourier 
transform (|7f|) the incident front can be considered as 
a superposition of numerous plane waves with energies 
larger and smaller than the height of the barrier whose 
interaction with the barrier can be examined indepen- 
dently due to above mentioned linearity of the problem. 
In spite of the fact that the mean energy of the front q is 
larger than the height of the barrier, the plane waves in 
the above superposition with energy less than the barrier 
height are completely reflected and this reflection results 
in the above mentioned overshoot. Next, we see the os- 
cillating after-effect in the local density of the reflected 
front (see the thick red curve in Fig. [3Ja)) that slowly 
tends to the pole contribution with increasing of £ shown 
by the thin black step-like curve. Comparing definitions 
2] and [3] and having in mind that both reflected wave 
components differ only by sign (v = —u) we conclude 
that the same oscillations are present in the current den- 
sity as well. This after-effect is closely related to the 
ZitterbewegungS. Usually this trembling motion of the 
electron is explained by means of interference between 
positive- and negative-energy relativistic wave function 
components. The interference itself is not sufficient and 
some additional perturbation is necessary. In our case 
this is the dispersive mass barrier that mixes the eigen- 
functions of the free electron motion. Compare with the 
reflection from the infinite barrier described in Sec. II 
where these components were not mixed. 



and calculated this integral numerically. A typical result 
is shown in Fig. [3Jb). 

Comparing the curves in both Figs. [3] (a) and (b) we 
note some differences. In the case of the below barrier 
reflection instead of an overshoot in the red thick solid 
curve that corresponds to the leading edge we see some 
diminishing of the electron density. This can be explained 
by the same superposition of many Fourier harmonics as 
in the previous above the barrier reflection case. Among 
them there are harmonics with energy larger than the 
barrier height. They do not reflect but penetrate the 
barrier what finally causes the above mentioned lack of 
density at the reflected front. 



V. PENETRATION INTO THE BARRIER 

The calculation of the front part that penetrates the 
barrier is similar to the calculation of the reflected one 
presented in the previous sections. The only difference 
is that the integral (|28p now is replaced by the integral 
(|29p. It leads to the different wave function components 
that depend on both (x and t) arguments. 

Calculating the integral (ffi?)) along the same contours 
in the complex fc-plane shown in Figs. [2ja) and (b) we 
divide this contribution into two parts. As in the case of 
reflection the pole contribution is trivial: 



e(t-x)T(q) I v?TT 



■2<1 



/ k I 



(34) 

It represents the result obtained by the scattering of 
plane waves truncated at the point x = t due to the 
finite velocity of the front. 

The contribution of the cut was calculated using the 
formulas analogous to Eqs. (|3"1?|) and (f3"3")) depending on 
whether the energy of the front is larger or smaller than 
the height of the barrier. 



B. Below the barrier reflection (q < 1) 

If the front energy is lower than the height of the bar- 
rier some problem with the cut contribution appears be- 
cause the pole is located exactly in the interval — 1 < x < 
1 and interferes with the integration contour Cb that is 
shown in Fig. [U(b). Therefore, we changed the cut and 
transformed the above integration contour into the Cd 
as shown in Fig. [U(c). The pole contribution gives the 
same result as in (f3"Tj) . Calculating the cut contribution 
as the integral along the Cd contour we changed the in- 
tegration variable k = exp[i(7r + tp)] in (|28p. presented 
this integral as 



V2 



tt/2 



dip^/cos tpc 



-£ cos tp 



-^(3^/2— £ sin tp) 

q + ie>-v 



-i(3<p/2— £ sin tp) 

q + ie-iv 



(33) 



A. Above the barrier penetration (q > 1) 

A typical result for the penetration of the front into 
the barrier with energy larger than the height of the bar- 
rier is shown in Fig. HJa). The local density and current 
is shown as a function of the coordinate x in the case of 
different time t values. As was already pointed out the 
barrier acts as a dispersive medium. That is why the 
form of the front isn't conserved, and the local density 
(actually the modulus of the wave function squared) and 
the current demonstrates a rather sophisticated behav- 
ior. Three important points are worth to be mentioned. 
First, the leading edge of the front is transformed into 
a sharp peak that moves with unit velocity inherent to 
the front in the dispersionless free motion region (x < 0). 
In the case of electromagnetic pulses it is known as a 
precursor—. The physical explanation of it's appearance 
is as follows. The theory of dispersion in the classical^ 
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Figure 4: (Color online) The local density Q (solid curves) 
and current (JS| (dotted curves) in the front that penetrates 
the barrier for different time t values. The thin vertical lines 
indicate the motion with the group velocity x = v £T t: (a) - 
above barrier reflection with q — 1.5 and (b) - below barrier 
reflection with q = 0.6. 



Second, as seen in Fig. IDJa) the main part of the front 
is retarded with respect to the above precursor. It caries 
the energy of the front which is the reason why it moves 
with the group velocity v gl - = de/dq = y^q 2 — 1/q. In the 
case of q = 1.5 we have v gv « 0.75, which is indicated by 
thin vertical lines in Fig.UJa). 

And at last the third point is that in spite of the ho- 
mogeneous barrier region we see the waving behavior far 
from the leading edge of the local density and current 
that are related to the above mentioned zitterbewegung, 
and which is much more pronounced in the penetrated 
front. 



B. Below the barrier penetration (q < 1) 

A typical result is shown in Fig. H{b) for the case of 
below barrier penetration when the front energy is lower 
than the height of the barrier. 

Here again we see the precursor moving with unit ve- 
locity that makes the below the barrier reflection anal- 
ogous to the above the barrier one. But now instead of 
preparing the barrier for the wave propagation this pre- 
cursor constructs step by step the evanescent wave that 
is inherent to the under the barrier reflection of plane 
waves. This process needs some resources. That is why 
the precursor loses its intensity (becomes narrower during 
its motion) in contrary to the penetration with the energy 
larger than the height of the barrier where such loses are 
not noticeable (compare the precursors in Figs, (a) and 
(b)). The oscillatory behavior of the penetrated fronts 
are also seen in the local density and current. This is in 
agreement with the fact that the mass barrier is a dis- 
persive medium. 



VI. PENETRATION OF THE RECTANGULAR 
WAVE PACKET INTO THE MASS BARRIER 



and quantum^ version, and the mechanical analogy of 
the Klein-Gordon equation^ as well, implies that the 
medium in which the waves propagate has its own degrees 
of freedom that actually causes the dispersion of waves in 
a stationary regime. This stationary regime corresponds 
to the excited state of the medium and, consequently, 
needs some time to be established. That is why the first 
piece of the front (the precursor itself) moves through the 
unexcited medium, and doesn't undergo any dispersion. 
Speaking figuratively one can imagine that the precursor 
prepares the media for the propagation of the main part 
of the pulse or front. In an analogous way one may con- 
sider the mass barrier region as a medium with inner de- 
grees of freedom whose quantum analog is the difference 
of the ground state energies of the two sublattices when 
we consider the graphene in the tight binding approxi- 
mation. Consequently, the precursor has to be present in 
our Dirac electron case. 



The reflection and penetration of the steep fronts con- 
sidered in the previous sections enables us to construct 
the result for the rectangular wave packet reflection by 
the mass barrier. The result is obtained as the super- 
position of two fronts with the proper shift A and am- 
plitudes chosen. The most interesting case is the below 
the barrier reflection that is shown in Fig. [5] in the case 
of q = 0.9. The reflected wave packet is shown in part 
(a). It depends on the argument £ = x + t as in the 
case of the front considered in the previous sections. The 
part of the wave packet that penetrates the barrier is a 
function of both x and t arguments. The local density 
is shown as a function of x in Fig. \5Hb) for given time 
t = 7 (red solid curve) together with the current (blue 
dashed curve). Here we see two precursors. One of them 
constructs the evanescent wave introducing some charge 
into the barrier, and later the other one destroys it ex- 
tracting that charge out of the barrier region. The part 
of that evanescent wave is clearly seen between those pre- 
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Figure 5: (Color online) (a) - the reflected rectangular wave 
packet with shift A = 5 between the opposite fronts as a 
function of £ = x + t: the real (blue dashed curve), imaginary 
(green dotted curve) wave function components and the local 
density (red solid curve); black dashed rectangular shows the 
wave packet reflected by a hard wall, (b) - the local density 
(solid curves) and current (dotted curves) for the penetrated 
part of the wave packet, at a given time. 



cursors. It is interesting to establish whether all charge is 
extracted or some of it is left in the barrier and continues 
its motion towards +oo. 

The easiest way to clear this point up is to compare 
the density in the reflected wave packet (the red solid 
curve in Fig. (Ha)) with black dashed curve in the same 
figure that represents the wave packet reflected by the 
hard wall. It is evident that the area below the red solid 
curve is smaller than the one below the black dashed 
curve. This means that the number of particles reflected 
by the finite mass barrier is smaller than this amount in 
the incident wave packet, consequently, some amount of 
the wave packet penetrates the barrier and moves there 
towards +oo in spite that its mean energy is smaller than 
the height of the barrier (q < 1). This fact is even bet- 
ter seen in Fig. [5] where the total number of particles in 



the reflected wave packet is shown as a function of time. 
Notice that N is smaller than A*o = 5 (the number of par- 
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Figure 6: (Color online) The number of particles in the re- 
flected wave packet as a function of time for different q val- 
ues in the case when this number for incident wave packet is 
iVo = 5. 

tides in the incident wave packet) in asymptotic region 
(for large time t), indicating explicitly that the closer q 
is to unity the larger the amount of the wave packet that 
penetrates the barrier. 

In Fig. OJa) we see one more interesting feature. It is 
a long waving tail that appears due to the above men- 
tioned degrees of freedom of the barrier. They are ex- 
cited by the incident wave packet and then radiate the 
energy even after the incident wave packet disappears. 
This is remarkable that mathematically this tail follows 
from the contribution of the cut to the reflected wave 
function . In the case of large £ values (what corre- 
sponds to the large time t asymptotic) the contribution 
to the integral is given by the contour edges and leads 
to the standard power type l/£ (non-exponential) be- 
havior, that appears in many systems characterized by a 
continuous spectrum, which for example were discovered 
in hydrodynamics^. 



VII. REFLECTION AND PENETRATION OF 
THE GAUSSIAN WAVE PACKET INTO THE 
MASS BARRIER 

Now let us consider the opposite case: the penetration 
into the mass barrier of a wave packet with a rather soft 
envelope. We consider the Gaussian wave packet replac- 
ing the initial condition by the following one: 

= _ ^„.,._ (35) 



Applying the same scheme as in section IIVI we present 
the above wave function as the Fourier integral ([55]) with 
the following Fourier transform: 



f(k) = (4Tra) 1 ^ C - a{q - k '> 2/2 . 



(36) 



Inserting this Fourier transform into Eq. (|25[) and replac- 
ing the exponent exp(ifcx) by the exponential solution of 
time dependent Dirac equation (|21I22[) as we did before 
we obtain the following integral representations of the 
incident Gaussian wave packet in the free motion region 
— oo < x < 



* (in) 0M) 



IT 



1/4 




2tt 3 /4 I 1 
its reflected part in the same region 

<7 1/4 / 1 \ 



(37) 



2tt 3 / 4 



-1 



(38) 



x / dkR(k)c-^ q - k ^/ 2 -' lk ^ x+t \ 
Jc 

and the one penetrated into the barrier (0 < x < oo) 



* (tr) (x,i) 



q 1 / 4 
: 2tt 3 / 4 



dke~ iR{k)x 



c 



I 1 + R{k) \ (9 _ fe) 2 /2+ifc(a ._ 4) 
\ 1 - R(k) ) 



(39) 



The calculation of these three integrals differs, how- 
ever, from the calculation of the previous ones in the 
case of the steep front or rectangular wave packet. The 
matter is that due to the exponent with momentum k 
squared the integrand has a more sophisticated singu- 
larity at infinity. That is why the integration contour 
can not be shifted to ±ioo, the method applied previ- 
ously fails, and we have to look for other possibilities to 
consider those integrals. We demonstrate the method of 
calculation in the case of the most simple integral (|3T|) 
for the incoming wave packet. Analyzing the argument 
of the single integrand exponent we see that the com- 
plex fc-plane can be divided into four sectors as shown 
in Fig. [7](a). In two of them, shown by green shadow- 
ing, the exponent increase for \k\ — > oo. Consequently, 
the integration contour C has to avoid these two sectors, 
and must be located completely in the two other white 
sectors. It is evident that going from — oo to oo this in- 
tegration contour has to cross the saddle point kg that is 
obtained by equating the derivative of the argument of 
the exponent to zero: 



Jt-[a(k-q) 2 /2-ik(x-t)) 



k—kt 



and 



= fco — q— i{x— t) = 0, 
(40) 



k = q + i(x - t). 



(41) 



(a) 






Figure 7: (Color online) (a) - Forbidden sectors in the com- 
plex fc-plane and the integration contour for the calculation of 
the integrals P7I39[) . (b) - Transformation of the integration 
contour in the case x < t. 



The value of the integral can be easily estimated by the 
saddle point method. We used the numerical integration 
over the contour Co shown in Fig. [7jb) by the horizon- 
tal red solid line. By the way the integral (|37| can be 
calculated analytically and it reads 



. , e -{*-t) /la 



( 7 ra) 1 /4 



*(*.*) = ^ 



what corresponds to the initial condition (|35|) that moves 
with unit velocity to the right conserving its form, as it 
should be in the dispersionless free motion region (— oo < 
x < 0). Thus in the case of the Gaussian wave packet the 
saddle point contribution to the integrals is the analog of 
the contribution of the pole in the case of a rectangular 
wave packet (or in the pure exponent case). 

Integrals (|3"8|) and (|3T)|) are more complicated due to 
the radical in the function R{k) in the integrands. Be- 
cause of these radicals there are branching points, and 
we have to make the cut in the complex fc-plane connect- 
ing them as it is shown in Fig. [TJb). Calculating these 
integrals it is important to check the relative position of 
the cut and the saddle point that according to (J4TJ) de- 
pends on the sign of [x — t). So, if x > t the horizontal 
contour Co crossing the saddle point lays above the cut, 
and consequently, only that contribution has to be taken 
into account. In the case of x < t as we see in Fig. |7[b) 
the contour Co is located below the cut. In this case the 
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contribution of the cut (namely, the integration over the 
contour Cb) has to be added. In the case of integrals 
(|38f and ([39| we performed the integration numerically 
taking into account the above mentioned features of the 
contours. 

A typical result for the above barrier reflection of the 
Gaussian wave packet is shown in Fig. [5] In part (a) 




Figure 8: (Color online) Above the barrier reflection of the 
Gaussian wave packet: the blue dashed and green dotted 
curves - the real and imaginary part of wave function compo- 
nent u (for reflected wave packet it coincides with the compo- 
nent —v), and the red solid curve - the density, (a) - reflected 
wave packet as a function of £ = x + t, the thin dashed black 
curve - the density of the incident wave packet when it is re- 
flected from the hard wall as in Sec. [II] (b) - The wave packet 
that penetrates the barrier, the thin dashed black curve - the 
density of the incident wave packet propagating in the absence 
of the barrier. 



crease) of intensity in the Gaussian wave packet is slow, 
the barrier manages to adjust itself for reflection, and the 
sophisticated features that we met in the case of rectan- 
gular wave packet reflection are not present. The single 
reminder of the previously considered reflection of the 
rectangular wave packet is the small asymmetry of the 
reflected Gaussian wave packet, some retardation of it as 
compared with the reflection of the Gaussian wave packet 
by the hard wall (shown by the thin dashed black curve), 
and the long waving tail after the reflected wave packet. 
The wave packet that penetrates the barrier is shown in 
Fig. [5](b) . This differs qualitatively from the rectangular 
wave packet case. There is no precursors, and the form of 
the penetrated wave packet is rather close to a Gaussian. 
It is remarkable that this form is more or less the same 
even in the case of below the barrier reflection (q < 1) 
as seen in Fig. [9ja) , only the amplitude being smaller 
and the width larger what is caused by the dispersion of 
the barrier medium. We see no precursors and no forma- 
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the reflected wave packet is shown. We see that it differs 
form the one of the rectangular wave packet or front (see 
for comparison Figs. [3] or [5{a)). There is no overshoot 
neither a lack of intensity at the edges, and the reflected 
wave packet conserves more or less its Gaussian form. 
The matter is that the Gaussian wave packet actually 
has no leading and trailing edges. The increase (or de- 



Figure 9: (Color online) (a) - the same as in[8][&) but for the 
below the barrier reflection case; (b) - the total number of 
particles in the barrier ([6]). 

tion of an evanescent wave. So, in the case of the below 
the barrier reflection the penetrated wave packet exhibits 
qualitatively the same Gaussian form with small asym- 
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metry and this packet is essentially retarded as compared 
with the motion of the wave packet in the absence of the 
barrier (shown by the thin dashed black curve), and a 
rather large spreading that, as it was already mentioned, 
appears due to the wave dispersion in the barrier. To 
our mind the retardation of the penetrated wave packet 
apparently illustrates the problem of the time interval 
that the wave packet spends in the barrier that was in- 
tensively discussed 15-20 years ago (see the discussion 
in Ref.[|l). 

In the case of the below the barrier penetration the 
part of the Gaussian wave packet continues its motion in 
the barrier towards +00 as it was in the rectangular wave 
packet case, what is evidently seen in Fig. EJb) where 
the total number of particles in the barrier is shown as a 
function of time t. This is clearly demonstrated by the 
non-zero asymptotic of these values at large t explicitly 
indicate that. These values increase when the energy q 
approaches the top of the barrier being essentially larger 
for short wave packets. 

VIII. CONCLUSIONS 

We considered a simple ID model of wave packet re- 
flection and penetration into a mass barrier. Interesting 
behavior was already found when considering reflection 
off a hard wall which is due to the fact that the wave 
packet isn't a point particle. For example we found al- 
ready reflection before the wave packet reached the wall, 
there was a compression of the wave packet during reflec- 
tion and the width of the wave packet was restored after 
reflection. 

The dependence on the wave packet shape was dis- 
cussed by considering two limiting shapes, namely, the 
rectangular wave packet with sharp leading and trailing 
edges, and the Gaussian wave packet which has a rather 
soft form. 

The most crucial differences of the wave packet reflec- 



tion and penetration into the mass barrier with those 
calculated for standard scattering of exponents were re- 
vealed in the case of the rectangular wave packet. Here 
the dispersion of the electron wave in the barrier showed 
itself to the full extent. We found different structures 
moving with different velocities. The leading edge of the 
front moves with velocity equal to unity because it moves 
through the unprepared dispersive medium of the bar- 
rier. Meanwhile the main piece of the wave packet (or 
front) moves, however, with the group velocity which is 
the reason why it is retarded. 

At the leading edge of the reflected front we see some 
overshoot or failure of some intensity due to the fact that 
the front itself is a superposition of many exponents that 
are reflected with different probabilities. In the case of a 
soft Gaussian wave packet the features of reflection and 
penetration are quite different. The main difference is 
the absence of the precursors. The matter is that the soft 
leading edge of the wave packet prepares the medium of 
the dispersive barrier gradually, and there is no possibil- 
ity for the precursor (moving in the unprepared media) 
to appear. Nevertheless some rudiments of the above fea- 
tures can still be revealed. That is the distortion of the 
wave packet form, making its leading edge sharper and 
the trailing edge more prolonged. Next, we see some re- 
tardation of the reflected and penetrated wave packet as 
compared with the motion in the absence of the barrier. 

We also demonstrated that for the Gaussian wave 
packet in the case when the mean wave packet energy 
is lower than the height of the barrier part of the wave 
packet penetrates into the barrier and moves along it, 
more or less conserving its shape. 
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